Dynamics of the self-interacting chameleon cosmology 
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Abstract 

In this article we study the properties of the flat FRW chameleon cosmology in which the cosmic 
expansion of the Universe is affected by the chameleon field and dark energy. In particular, we 
perform a detailed examination of the model in the light of numerical analysis. The results illustrate 
that the interacting chameleon filed plays an important role in late time universe acceleration and 
phantom crossing. 
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1. INTRODUCTION 



Recently, the observations of high redshift type la supernovae and the surveys of clusters 
of galaxies jl]-^] reveal the universe accelerating expansion and that the density of matter 
is very much less than the critical density. Also the observations of Cosmic Microwave 
Background (CMB) anisotropies indicate that the universe is flat and the total energy density 
is very close to the critical one j^J. 

The above observational data properly complete each other and point out that the dark 
energy (DE) is the dominant component of the present universe which occupies about %73 
of the energy of our universe, while dark matter (DM) occupies %23, and the usual baryonic 
matter about %4. There are prominent candidates for DE such as the cosmological con- 
stant , a dynamically evolving scalar field ( like quintessence) [8j, |9| or phantom (field 
with negative energy) 10[ that explain the cosmic accelerating expansion. Meanwhile, the 



accelerating expansion o: 
cosmology and so on [12 



universe can also be obtained through modified gravity 



11], brane 



Two of the most serious issues with regards to the DE models, in particular with cosmo- 
logical constant as a candidate, are the fine tuning problem and cosmic coincidence problem. 
The absence of a fundamental mechanism which sets the cosmological constant to zero or 
very small value is the cosmological constant "fine-tuning" problem and a good model should 
limit the fine tuning as much as possible. The problem of comparability of the DE density 
and the DM energy density at the recent epoch known as the coincidence problem and one 
of the most frequently used approach to moderate the cosmological coincidence problem is 
the tracker field DE scenario [9|. The DE can track the evolution of the background matter 
in the early stage, and only recently, it has negative pressure, and becomes dominant . Thus, 
its current condition is nearly independent of the initial conditions 19]- 23]. 

On the other hand, to explain the early and late time acceleration of the universe, it is 
most often the case that such fields interact with matter; directly due to a matter Lagrangian 
coupling, indirectly through a coupling to the Ricci scalar or as the result of quantum loop 



corrections 



24j- 28] . If the scalar field self-interactions are negligible, then the experimental 



bounds on such a field are very strong; requiring it to either couple to matter much more 



weakly than gravity does, or to be very heavy 



29]- 32]. Unfortunately, such scalar field is 



usually very light and its coupling to matter should be tuned to extremely to small values 



2 



in order not to be conflict with the Equivalence Principal 33]. 

An alternative attempt to overcome the problem with light scalar fields has been suggested 
in chameleon cosmology 34[-36(. In the proposed model, a scalar field couples to matter 
with gravitational strength, in harmony with general expectations from string theory whilst 
at the same time remaining very light on cosmological scales. The scalar field which is 
very light on cosmological scales is permitted to couple to matter much more strongly than 
gravity does, and yet still satisfies the current experimental and observational constraints. 
The cosmological value of such a field evolves over Hubble time-scales and could potentially 
cause the late-time acceleration of our Universe {stJ. The crucial feature that these models 
possess are that the mass of the scalar field depends on the local background matter density. 



While the idea of a density- dependent mass term is not new 
in 



38 



43J, the work presented 



36] 371 ] is novel in that the scalar field can couple directly to matter with gravitational 



strength. 



Considering the possible interaction between DE and background matter [44|], the whole 
system (including the background matter and DE) may be eventually attracted into the 
scaling attractor, a balance achieved, thanks to the interaction. In the scaling attractor, 
the effective densities of DE and background matter decrease in the same manner with the 
expansion of our universe, and the ratio of DE and background matter becomes a constant. 
So, it is not strange that we are living in an epoch when the densities of DE and DM 
are comparable. In this sense, the cosmological coincidence problem is moderated. The 
dynamical attractor of the cosmological system has been employed to make the late time 
behaviors of the model insensitive to the initial condition of the field and thus moderates 
the fine tuning problem. 

In this paper, we assign two important roles to the chameleon scalar field to describe the 
late time acceleration of the universe and possibly predict the fate of the universe. Its first 
role as already expressed is that the chameleon field which is very light on cosmological scales 
and its mass strongly depend on the local background matter density satisfies the current 
experimental and observational constraints. Its second role is that since the field mimics the 
background (radiation/matter) matter field, subdominant for most of the evolution history 
except at late times when it becomes dominant, it may be regarded as a cosmological tracker 
field. Tracker models are independent of initial conditions used for field evolution but do 
require the tuning of the slope of the scalar field potential. During the scaling regime the 



3 



field energy density is of the same order of magnitude as the background energy density. 



This work is different from that of ref. 



45] in which we assume that the coupling function 



to the matter is exponential while in there a linear coupling to the matter is assumed. Our 



work also differs from that of ref. 



37] in that the potential function for the chameleon field 



in there is in exponential form that in its series expansion is a constant plus corrections in 
terms of a stability parameter. 



2. THE MODEL 



In this section we consider the chameleon gravity with the action, 

S = f - \<t>,^ + V(<j>) + f(<f>)L m }^dx\ (1) 

where R is Ricci scalar, G is the newtonian constant gravity and <fi is the chameleon scalar 
field with a potential V{4>). The modified matter is f(<p)L m , where f{<p) is an analytic 
function of <fi and describe the nonminimal interaction between the matter and chameleon 
field. In a spatially flat FRW cosmology, the variation of action ([T]) with respect to the 
metric tensor components leads to the field equations, 

3H 2 = p m f+±<j) 2 + V (</>), (2) 
2H + 3H 2 = - 1 -<P 2 + V(<P), (3) 

where we put 8irG = 1 and H = - with a the scale factor. We assume a perfect fluid 
with p m = 7p m where p m is the contribution from the matter to the energy density. In 
the following we assumed that the matter in the universe is cold dark matter (CDM) with 
7 = 0. Variation of the action ([Q) with respect to scalar field gives the wave equation for 
chameleon field as, 

4> + 3H<P=-V' -\ Pm f\ (4) 

where prime indicated differentiation with respect to (p. From equations (j2J), (J3J) and 01]), 
one can easily arrive at the relation, 

(Pmf) + 3Hp m f = -pm&f'. (5) 
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From equations fl2]) and ([3]) and in comparison with the standard friedmann equations we 
identify p e ff and p e ff as, 

Peff = Pmf + t^ 2 + V(4>) = p ch + p de , (6) 

Peff = \ft - V {4>) = Pde, (7) 

where "de" and "ch" stand for dark energy and chameleon, respectively and p e ff = u e ffp e ff. 
In here, we introduced p c h = p m f and call it "chameleon energy density". The conservation 
equation for dark energy and chameleon field coupled to matter separately are, 

p ch + 3Hp ch = Q, (8) 
p de + 3H(l +u d )p de = -Q, (9) 

where Q is the interaction term and p de = u de p de . Comparing equations ([5]) and ([8]), we find 
that Q = \pm<\>f' ■ By defining jM = g(t, (p, 0) we can rewrite Q as 

Q = \gPeh = \g^\~l + We//) - 2V(1 + UJeff))- (10) 

From equation ffTUl) . one observe that at t — t cross when u e ff = — 1, then 

Q = -\g<P 2 - (ii) 

In Q, the variable g gauges the intensity of the coupling between matter and chameleon 
field. For g = 0, there is no interaction between chameleonic dark matter and dark energy. 
The Q term measures the different evolution of the DM due to its interaction with DE which 
gives rise to a different universe expansion. From equation (jSJ) and ffTU]) . derivative of Q with 
respect to N = ln(a) gives, 

with the solution, 

Q = Q ga~ 3 ey^ dt = Q ga' 3 f^. (13) 



One can also rewrite equation ([12]) in terms of redshift as 
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By defining the ratio of the chameleon energy density to dark energy density as r = — , we 

Pde 



then obtain, 



1 3H 2 

1 + - = . 15 

r pch 



From equations (jl5l) and ([8]) one finds 



r 3rH 2 #, Q 

+ ^-:( 1 +o77o) = T^- (16) 



(1 + r) 2 1 + r 3H 2 ' 3H 2 ' 

In addition, using u e ff = — 1 — |-^-, we obtain 

r 3rHu eff Q 



(1 + r) 2 1 + r 3H 2 ' 
Using z — — (1 + z)H(z), one can rewrite equation ( 1T71) in terms of the redshift z as, 

(18) 



'1 + z) dr 3rco ef f _ Q 



(l + r) 2 dz 1 + r 3H 3 
A numerical discussion of the model is presented in the next section. 



3. NUMERICAL DISCUSSION 



A numerical analysis is performed in the following to discuss the model in terms of the 
dynamical variables. Using equation ( TTT]) . Fig.l, shows the dynamic of the EoS parameter, 
u) e ff, the interaction term Q and —^gcj) 2 . As illustrated, the phantom crossing occurs at 
two different situations, for Q ^ at t ^ when Q = —\g4> 2 ^0. Or alternatively for 
Q = —\g4> 2 = 0, the phantom crossing does not occur. Thus, phantom crossing may occur 
only when there is an interaction between chameleon and dark energy. 



To better understand the cause of phantom crossing i the model we consider two cases 
in the interaction between the chameleon and dark energy fields: 1) the interaction term 
is Q = —\g<fi 2 with exponential coupling function as /(</>) = foe x ^, 2) the interaction term 
arises from holographic dark energy. 



case 1 

If we assume that f(<p) = foe x< ^ where A is constant, then 



9 = A0, (19) 
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Fig.l: The dynamic of effective EoS parameter, uj e ff, the interaction term Q and 
and from equation (TIDl) . we rewrite Q as 



1 „J,2 



Q = + We//) " + ^ )} - (20) 



Fig.2 shows the dynamics of the effective EoS parameter, u) e ff, interaction term Q and 

a 3 



3 

— jW for different values of A. As expected, the crossing happens only when non vanishing 



3 

Q and —j^ intersect each other. For A = —2 the phantom crossing occurs at t — —0.355 
and t = 0.45. The graphs show that for A decreasing to —3 and —5, the gap between two 
phantom crossing location reduces. Eventually, for A — > — oo or Q&zg —> 0, the two points 
merge into one and the effective EoS parameter becomes tangent to the divide line (see 
Fig.3). 



case2 

If we assume that the interaction term is give by Q = 3aHp c h where a is interacting 
parameter, then by comparing it with the interaction term in ( |T0|) . we yield 

9 = j = l^aH (21) 

where a = Furthermore, from equation([8]) we can obtain chameleon energy density as 

Pch = Poa- 3{1 - a) = p a- 3 a 3a , (22) 

where p c h increases by a factor of a 3cr . Obviously, in case of no interaction between dark 
sectors we have p c h = Po a ~ 3 - Solving equation (121!) for /, we find / = foa 12cr where for 
different values of a, we can obtain the dynamic of the reconstructed function /. Using the 
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Fig. 2: The evolution of effective EoS parameter, w e //, 

3 

interaction term Q and — |w for A = —2, —3, —5. 



Fig.3: The evolution of effective EoS parameter, w e //i for A — > — oo or Qhg — > 



above equations, after some calculation, we obtain 

1 dg 3 

or in terms of redshif z 



(l + w e// ) 



-{L + U e ff) 



(23) 



(24) 



g dz 2 

Fig. 4 shows the dynamic of the effective EoS parameter in comparison with From the 
graph we see that phantom crossing occurs when ^ vanishes. 



S 



dg(z) 



Fig.4: The evolution of effective EoS parameter w e // 5 and 



Using equation (I24p . equations (TT4")) for Q and (121 p for redshift 2; can be rewritten as 



^ = Q[ g \ 9 1 3 ggg ] 
cfe " V1 4i 2(1 + 2) 2(1 + z) ] 

= g(z)(l + z) 

12a 



(25) 
(26) 



Inserting fl26|) into fl25|) . we obtain 

rfQ 3g 11 

^ = TT7 ( - (T+ 2 + 2^ ) - (27) 

Using ffTB"]) . from equation (1271) we find a dynamical equation for the ratio of the chameleon 
energy density to dark energy density in terms of interaction parameter o as 



1 + z dr 
r(l + r) dz 

Or, using equation (l2~Tj) we have, 

dQ „ r 9(l-(7 



3(a + w e// ) 



(28) 



} 



(29) 



dz ^2 1 + 2 2r(r + l)rf2 
Using equation ({TBI) we can solve the differential equation (129]) for Q to find the solution: 

g = C(^)^(i + 2)K 1 -^ (30) 

where C is a constant. Also we also drive the parameter g in terms of ratio r and as the 
interaction parameter a 



C r + l 1 
= t( W 1 



(31) 



In Fig. 5 we depicted the dynamics of the effective EoS parameter for positive and negative 
values of a. 
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Fig. 5: The evolution of effective EoS parameter 0J e ff f° r different a 



For positive values of a, it shows that at z = 0.43 in the past the effective EoS parameter 
vanishes. The phantom crossing occurs both in the past and future and the universe enters 
accelerating phase in the near past. For some negative interacting parameters, for example 
a = —1/2, the effective EoS parameter is always negative in the past which is not supported 
by the observational data. However, interestingly, for a = — 1, the graph shows that the 
universe undergoes acceleration at about z ~ 0.4 in near past whereas phantom crossing 
occurs in near future. In addition the current value of EoS parameter is about ou e ff ~ —0.7 
and also approaches zero in the far past which is strongly favored by the observational data. 



4. SUMMARY 



In this work, we study the dynamics of the universe in chameleon cosmology. In detailed 
examination, we investigate the interacting phantom cosmological paradigm where the scalar 
field acts as dark energy and interact with chameleon field taken as the coupling /(</>) to 
matter lagrangian. By introducing the function g(t, <f), <fi) as the relative ration in the function 
/(</>), we express all the dynamical variables such as the interacting Q, the ratio of the 
chameleon energy density to dark energy density, r and the effective EoS parameter u e ff 
in terms of or in relation to g. In two distinguished scenarios, with numerical analysis we 
illustrate the conditions for universe acceleration and phantom crossing. In particular, in 
the second scenario, for some specific values of interacting parameter er, the model predicts 
universe acceleration and phantom crossing in the near past, in addition to recovering matter 
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dominated universe in the far past which is supported by observational data. 
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